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EXECUTIVE  SUMMARY 


PROBLEM 

A  method  is  needed  for  finding  the  variability  of  the  electric  field  strength  in  the 
earth-ionosphere  waveguide  at  very  low  frequencies  due  to  variations  in  the  electron 
density  profile. 

RESULTS 

A  method  is  developed  for  finding  the  standard  deviation  of  the  field  strength  in  the 
earth-ionosphere  waveguide.  This  method  uses  derivatives  of  the  waveguide  eigenangles 
with  respect  to  height  and  slope  variations  in  the  electron  density  profile.  The  calculated 
results  are  in  reasonable  agreement  with  available  data. 
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1.  Comparison  of  calculated  values  with  GBR  data.  Frequency  =  16.0 
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2.  Comparison  of  calculated  values  with  Anthom  data.  Frequency  = 
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field  strength  =  0.48  gauss.  Parameters  of  the  profile  are 
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h'  =  85.3  km  and  =  0.5  km'* .  12 
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INTRODUCTION 


Predictions  of  the  field  strengths  of  signals  from  very  low  frequency  (VLF)  and  low 
frequency  (UF)  transmitters  are  of  considerable  interest  to  the  Navy.  For  example,  the 
choice  of  frequency  diversity  depends  on  the  locations  of  nulls  in  the  signal  level.  The 
predictions  of  field  strengths  are  based  on  mode  theory  for  the  earth-ionosphere 
waveguide  [Budden,  1961;  Pappert,  Gossard,  and  Rothmuller,  1967;  Morfitt  and 
Shellman,  1976].  Variations  occur  in  the  electron  density  profile,  however,  and  a  method 
of  routinely  calculating  the  standard  deviation  of  the  signal  strength  as  a  function  of 
distance,  due  to  these  variations,  has  not  been  available.  Such  a  method  is  needed  for 
more  realistic  assessments  of  signal  coverage. 

In  this  report  an  exponential  form  of  the  electron  density  is  assumed.  The  overall 
height  is  given  by  h'  and  the  slope  by  ^  such  that 

N  =  1.43  X  10’  exp{/?(A-A')  -  0.15A] 

where  height  is  in  kilometers  and  electron  density,  V,  is  in  electrons  per  cubic  centimeter. 
The  parameters  h'  and  are  assumed  to  have  Gaussian  distributions  with  standard 
deviations  Oh  and  op  with  respect  to  time  variations  in  the  ionosphere. 

Ferguson,  Morfitt,  and  Hansen  [1985]  used  measurements  of  signal  levels  of  transmis¬ 
sions  from  NGR  in  Greece  at  59  kHz  to  deduce  h'  and  fi,  their  standard  deviations,  a* 
and  Op,  and  the  correlation  between  them.  The  measurements  were  made  at  two  locations 
on  either  side  of  a  modal  interference  null  at  900  km  and  1300  km  west  of  the 
transmitter.  The  object  of  this  present  work  is  to  derive  an  efficient  way  of  calculating  the 
standard  deviation  of  the  field  strength  as  a  function  of  distance  for  any  given  values  of 
Oh  and  Op.  It  is  assumed  that  there  is  no  correlation  between  h'  and  jS.  This  latter 
assumption  is  not  fundamental  to  the  formulation,  however,  which  could  be  expanded  to 
include  correlation. 

The  approach  taken  in  this  present  work  is  to  find  the  rate  of  change  of  the  sines  of 
the  waveguide  eigenangles  and  from  that  information  to  find  the  expectation  values  of  the 
field  strength  and  the  standard  deviation.  The  excitation  factors  are  assumed  not  to  vary. 

THEORY 


The  expectation  of  the  variance  is 

<  (£  -  <E>)  {E*  -  <£•>)> 

=  <££*>  -  <£<£•»  -  «£>£*>  +  <£><£*> 
=  <££•>  -  <£><£•> 
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where  the  electric  field,  £,  as  a  function  of  h'  and  0,  is  taken  to  be 


^  +  (dS:,/dh’)M'  +  (dS„/d0)A0]x} 

where  A„  is  the  excitation  factor,  S„  is  the  sine  of  the  earth-ionosphere  eigenangie,  and  k 
is  the  wave  number.  The  variation  with  respect  to  h'  and  0  is  taken  to  be  independent. 

In  general 

</(w)  >  =  J/(w)  pi(w)  dw 
where  the  integral  is  taken  from  -«  to  «  and 

Pi{w)  =  [l/a(2?r)*/2]  exp  (-m^/2o^) 


For  a  sum 


For  a„  of  the  form 


X  J  a«(w)pi(w)dw  =  X  <  > 


a„  =  exp  (-  Qo  -  Oiw) 

the  expectation  is 

<a„>  =  {\/a{2nyf'^]  j  exp(-ao  -  OiW  -  dw 

where  the  subscript  n  has  been  omitted  from  the  constants  oo  and  Ui  for  simplicity. 
Completing  the  square 

Qiw  +  w^/2o^  =  (cio^  +  wyii&  -  d\(yil 


so  that 

<a„>  =  [l/(7(27r)*/^]  exp(-ao  +  aiO^/2)  |  exp(-(aiO^  +  w)V2o^]dw 
The  integral  is  equal  to  so  that 

<a„>  =  exp(-  oo  +  0 10^/2) 

The  result  is  easily  extended  to  the  case  of  two  independent  variables,  Wa  and  w^,;  that 
is,  for 


a„  =  exp(-ao  -  oiw^  -  b\Wbi 


2 


The  expectation  is  then  of  the  form 


where 


<a„> 


Wb)p2(Wa  ,  Wb)dWa  dWt 


P2{Wa  .  Wb)  =  (l/2jr  a^Ob)  exp[-  (wi/o^  +  wl/ai)/2] 

so  that 

* 

<  a„  >  =  exp(-  oo  +  ala^/l  +  h\db/l) 

The  expectation  of  £  is  then  proportional  to 

<£>  —  ^Xn  exp  [-ik  S,^-(kxak  dSJdh'Y/2  -  {kxa^  dSJd^^/2\ 

Note  that  the  expectation  of  £  is  not  the  value  of  £  for  AA'  =  =  0.  This  is  because 

of  the  asymmetrical  nature  of  the  exponential  function. 

The  expectation  of  EE*  is  proportional  to 

<EE*>-  exp{-/^(S„  -  Sl,)x-[kxab{bSJdW  -  dSl/dh’)Y/2 

-  [kxap{dS„}dfi  -  dS;/d^]V2> 

The  derivatives  of  sine  0  with  respect  to  h'  and  ^  were  computed  by  finding 
eigenangles  for  finite  increments  in  h'  and  ^ .  The  independence  of  h'  and  /8  is  assumed. 
The  derivatives  are  then 

dSi/bh'  =  cos6>{[0,(A'  +  -  AA')]/2AA'} 

=  cos0([0,O5  +  A/5)  -  Oiifi  -  A^]/2A^ 

Only  the  10  or  12  most  significant  modes  were  included.  Care  was  taken  to  ensure  that 
the  same  set  of  modes  was  used  for  the  increments  and  decrements  in  h'  and  ^  as  for  the 
^  no-increment  set.  Increments  isJi'  =  1km  and  A/?  =  0.1  were  used  for  this  purpose. 

In  principle,  analytic  derivatives  of  the  eigenangles  could  be  computed  so  that  ordering 
#  would  not  be  a  problem.  Although  computer  code  for  this  is  already  available,  it  is  large 

and  would  unduly  complicate  the  MODESRCH  program. 

RESULTS 

Calculations  of  the  variability  of  electric  field  strength  in  the  earth-ionosphere 
waveguide  were  compared  with  data  taken  on  board  a  ship  in  the  Atlantic,  the 
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GTS  Callaghan,  from  which  GBR  at  Rugby,  England;  the  Anthom  transmitter  in  England; 
NSS  at  Annapolis;  and  NAA  at  Cutler,  Maine,  were  monitored  [Computer  Sciences 
Corporation,  1989].  The  data  were  taken  at  night  on  many  traverses  of  the  Atlantic,  and 
hence  the  superimposition  of  these  data  gives  a  good  indication  of  the  nighttime 
variability  in  the  £-field  due  to  variation  in  the  ionosphere. 

In  the  calculations,  each  path  was  approximated  to  be  entirely  over  seawater  even 
though  in  fact  a  small  part  of  each  path  in  the  vicinity  of  the  transmitter  was  over  land. 
Each  path  was  approximated  to  be  homogeneous  in  other  respects  as  well.  Actually  the 
azimuth  of  propagation  and  the  dip  and  strength  of  the  earth’s  magnetic  field  varied 
somewhat  along  each  path.  For  seawater,  the  relative  permitivity  €r  =  81,  and  the 
conductivity  ct  =  4.0.  Twelve  modes  were  used  for  comparisons  with  NSS  data  and  ten 
modes  were  used  for  each  of  the  other  cases. 

Figures  la,  2a,  3a,  and  4a  show  a  comparison  between  field  strength  data  and  field 
strength  calculations  for  GBR  (16  kHz),  Anthom  (19  kHz),  NSS  (21.4  kHz),  and  NAA 
(24  khz),  respectively  (figures  are  at  the  end  of  the  report).  The  fit  is  fairly  close, 
considering  that  homogeneity  along  the  path  is  assumed. 

Figures  lb,  2b,  3b,  and  4b  show  the  data  with  the  expectation  of  the  field  strength  as  a 
function  of  distance.  Curves  of  standard  deviation  are  also  shown.  These  are  for  a*  =  1.0 
and  ap  =  0.05.  It  is  seen  that  the  magnitudes  of  the  standard  deviations  for  the  computed 
curves  are  about  the  same  as  for  the  data. 

The  expectation  of  the  field  strength  and  the  standard  deviations  are  also  shown  in  the 
remaining  figures  for  various  combinations  of  a*  and  ojj.  In  figures  Ic,  2c,  3c,  and  4c, 
Oh  =  1.0  and  ojj  =  0.05,  and  the  usual  curve  of  field  strength  versus  distance  is  shown  as 
a  dashed  line. 

In  figures  Id,  2d,  3d,  and  4d,  ~  1.0  and  ojj  =  0.0,  and  in  figures  le,  2e,  3e,  and 
4e,  a*  =  0.0  and  a^s  =  0.05.  These  last  figures  show  the  separate  effects  of  variation  in 
h'  and  The  variation  of  ^  affects  the  field  mostly  at  long  distances.  This  would 
correspond  to  the  fact  that  the  slope  of  the  electron  density  profile  affects  the  attenuation 
much  more  than  does  the  height. 

CONCLUSIONS 

A  method  has  been  presented  of  finding  the  standard  deviation  of  the  field  strength  in 
the  earth-ionosphere  waveguide.  This  method  uses  derivatives  of  the  eigenangles  with 
respect  to  height  and  slope  variations  in  the  electron  density  profile.  The  calculated  results 
are  in  reasonable  agreement  with  available  data. 

A  method  that  uses  the  derivatives  of  the  eigenangles  is  considered  to  be  better  than 
one  that  uses  the  derivatives  of  E  or  ln|£|  directly  since  the  variation  of  the  eigenangles  is 
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more  nearly  linear.  The  curve  of  the  expectation  of  ln(£|  is  smoother  than  that  of  lnl£(, 
and  the  expectation  of  the  standard  deviation  of  ln|£(  shows,  not  surprisingly,  a 
considerable  amount  of  asymmetry. 

The  method  has  some  drawbacks.  Sets  of  modes  for  the  various  profiles  needed  to 
find  the  derivatives  must  be  visually  checked  for  each  case,  and  the  situation  is 
ambiguous  when  modes  are  close  together.  Variation  of  the  excitation  factors  is  not 
accounted  for.  Furthermore,  the  amount  of  computation  is  proportional  to  the  square  of 
the  number  of  modes  once  the  sets  of  modes  are  found.  Nevertheless,  the  method  is 
expected  to  be  useful  in  predicting  the  variation  of  signal  strength  in  the  earth-ionosphere 
waveguide. 
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(a)  Calculated  field  strength  compared  with  GBR  data. 


(b)  Expectation  of  the  field  strength  and  standard  deviation 
compared  with  GBR  data,  a*  *  1.0  km  and  ap  »  0.05  km‘^  . 

Figure  1.  Comparison  of  calculated  values  with  GBR  data. 
Frequency  *  16.0  kHz,  azimuth  =  265  degrees,  dip  » 

66  degrees,  and  magnetic  field  strength  =  0.47  gauss. 
Parameters  of  the  profile  are  h'  =  83.8  km  and 
P  *  0.5  km"’  . 
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(c)  Field  strength  (dashed  curve),  expectation  of  field 
strength,  and  standard  deviation,  o*  =  1.0  km  and 
ap  -  0.05  km“^  . 


0.00  1.00  2.00  3.00  4.00  5.00  6.00  7.00  8.00 


Distance  (Mm)  (Contd) 

(d)  Expectation  of  the  field  strength  and  standard  deviation 
for  variability  of  the  height  of  the  profile.  Oh  =  1.0  km  and 

ap  =  0.0. 

Figure  1.  Continued. 


7 


ein 


140 


120  -I 


0.00  1.00  2.00  3.00  4,00  5.00  6.00  7.00  8.00 


Distance  (Mm) 

(e)  Expectation  of  the  field  strength  and  standard  deviation 
for  variability  of  the  slope  of  the  profile,  a*  *  0.0  and 
=  0.05  kin”’  . 

Figure  1.  Continued. 
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(a)  Comparison  of  calculated  values  with  Anthom  data. 


(b)  Expectation  of  the  field  strength  and  standard 
deviation  compared  with  Anthom  data,  =  1.0  km 
and  ap  =  0.05  km"’  . 


Figure  2.  Comparison  of  calculated  values  with  AnUiom 
data.  Frequency  =  19.0  kHz,  azimuth  «  265  degrees,  dip 
=  69  degrees,  and  magnetic  field  strength  *  0.48  gauss. 
Parameters  of  the  profile  are  h'  »  84.0  km  and 
P  =  0.5  km"’  . 
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(c)  Field  strength  (dashed  curve),  expectation  of  field 
strength,  and  standard  deviation,  a*  «  1.0  km  and 
ojj  =  0.05  km“*  . 


(d)  Expectation  of  the  field  strength  and  standard  deviation 
for  variability  of  the  height  of  the  profile,  a*  *  1.0  km  and 
ojj  «  0.0. 

Figure  2.  Continued. 
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E  in  nV/m  (dB) 


(e)  Expectation  of  the  field  strength  and  standard  deviation 
for  variability  of  the  slope  of  the  profile,  a/,  =  0.0  and 
Ofi  =  0.05  km“*  . 


Figure  2.  Continued. 


(a)  Comparison  of  calculated  values  with  NSS  data. 


(b)  Expectation  of  the  field  strength  and  standard 
deviation  compared  with  NSS  data,  oh  =  1.0  km 
and  Ofi  =  0.05  km'‘  . 


Figure  3.  Comparison  of  calculated  values  with  NSS 
data.  Frequency  =  21.4  kHz,  azimuth  =  85  degrees,  dip 
=  72  degrees,  and  magnetic  field  strength  =  0.53  gauss. 
Parameters  of  the  profile  are  h’  =  85.3  km  and 
P  =  0.5  km-‘  . 
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140  -j 
120  ■ 


Distance  (Mm) 

(c)  Field  strength  (dashed  curve),  expectation  of  field 
strength,  and  standard  deviation,  oi,  =  1.0  km  and 
ap  =  0.05  km"^  . 


(d)  Expectation  of  the  field  strength  and  standard  deviation 
for  variability  of  the  height  of  the  profile,  o*  *  1.0  km  and 
ap  =  0.0. 

Figure  3.  Continued. 
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E  jn  nV/m  (dB) 


(e)  Expectation  of  the  field  strength  and  standard  deviation 
for  variability  of  the  slope  of  the  profile,  as  ~  0.0  and 
ap  -  0.05  km'^  . 

Figure  3.  Continued. 
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(a)  Comparison  of  calculated  values  with  NAA  data. 


(b)  Expectation  of  the  field  strength  and  standard 
deviation  compared  with  NAA  data,  a*  =  1 .0  km 
and  -  0.05  km“^  . 


Figure  4.  Comparison  of  calculated  values  with  NAA 
data.  Frequency  =  24.0  kHz,  azimuth  =  120  degrees,  dip 
®  70  degrees,  and  magnetic  field  strength  =  0.50  gauss. 
Parameters  of  the  profile  are  h'  =  85.1  km  and 
0  =  0.43  km-‘  . 
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0.00  1.00  2.00  3.00  4.x  5.X  6.X  7.X  8.X 

Distance  (Mm) 

(c)  Field  strength  (dashed  curve),  expectation  of  field 
strength,  and  standard  deviation,  =  1.0  km  and 
Ofi  =  0.05  km"*  . 


(d)  Expectation  of  the  field  strength  and  standard  deviation 
for  variability  of  the  height  of  the  profile,  o),  *  1.0  km  and 

Op  =  0.0. 


Figure  4.  Continued. 
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(e)  Expectation  of  the  field  strength  and  standard  deviation 
for  variability  of  the  slope  of  the  profile,  a*  =  0.0  and 
Ofi  =  0.05  km'^  . 

Figure  4.  Continued. 
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